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Abstract 


The electromagnetic scattering from an arbitrarily shaped aper- 
ture backed by a rectangular cavity recessed in an infinite ground 
plane is analyzed by the integral equation approach. In this approach , 
the problem is split into two parts: exterior and interior. The electro- 
magnetic fields in the exterior part are obtained from an equivalent 
magnetic surface current density assumed to be flowing over the aper- 
ture and backed by an infinite ground plane. The electromagnetic 
fields in the interior part are obtained in terms of rectangular cavity 
modal expansion functions. The modal amplitudes of cavity modes are 
determined by enforcing the continuity of the electric field across the 
aperture. The integral equation with the aperture magnetic current 
density as an unknown is obtained by enforcing the continuity of mag- 
netic fields across the aperture. The integral equation is then solved 
for the magnetic current density by the method of moments. The elec- 
tromagnetic scattering properties of an aperture backed by a rectangu- 
lar cavity are determined from the magnetic current density. 

Numerical results on the backscatter radar cross-section (RCS) pat- 
terns of rectangular apertures backed by rectangular cavities are com- 
pared with earlier published results. Also numerical results on the 
backscatter RCS patterns of a circular aperture backed by a rectangu- 
lar cavity are presented. 

1. Introduction 

Electromagnetic (EM) scattering characteristics of metallic cavities are useful in studying radar 
cross section and electromagnetic penetration properties of objects consisting of these cavities as sub- 
structures. A large body of work which treats these cavities has been reported in the literature (e.g., 
refs. 1 to 12). In references 9 and 10, the problem of determining the RCS of an aperture backed by a 
cavity recessed in an infinite ground plane has been solved by the equivalence principle and integral 
equation approach. The integral equation in references 9 and 10 was solved by discretizing the aperture 
into rectangular or quadrilateral patches for use in the method of moments. However, when the aperture 
is of irregular shape with sharp comers and bends or loaded with microstrip patches of irregular shape, 
the discretization scheme used in references 9 and 10 may not be able to model the aperture efficiently. 
Discretization of an irregular-shaped aperture into triangular subdomains may be preferred over the 
quadrilateral patches. 

The problem of an arbitrarily shaped aperture backed by a rectangular cavity is formulated in terms 
of an integral equation with the magnetic surface current density in the aperture as an unknown. Dis- 
cretizing the aperture into triangular subdomains and using the method of moments transformed the 
integral equation into a matrix equation, which is solved by standard lower or upper (LU) decomposi- 
tion techniques. From the magnetic surface current density, the scattered far fields from the arbitrarily 
shaped aperture are determined. The present method is validated by comparing the numerical results 
obtained by the present formulation with the results available in reference 9. 

The formulation of the problem in terms of integral equations using the surface equivalence princi- 
ple is developed in section 2, numerical results on the RCS of rectangular apertures backed by rectangu- 
lar cavities are presented and compared with earlier published results in section 3, and the RCS of a 
circular aperture backed by a rectangular cavity in an infinite ground plane is also presented in 
section 3. 



2. Theory 


Consider an aperture backed by a rectangular cavity recessed in an infinite ground plane as shown 
in figure 1(a). The cavity is assumed to be illuminated by an electromagnetic plane wave. The problem 
may then be divided into interior and exterior parts as shown in figures 1(b) and 1(c). With the equiva- 
lence principle, the fields in the exterior region can be found by replacing the entire ground plane 
including the aperture with equivalent electric surface current density j and the aperture with equivalent 
magnetic surface current m only. However, with the application of image theory, only the magnetic cur- 
rents are needed to determine the exterior fields. From the boundary condition of continuous electric 
fields, the tangential electric fields inside are expressed in terms of-m. If E w and H w are the EM fields 
due to -m inside the cavity and E s , and H v are the EM fields due to m in the exterior region, then using 
the continuity of tangential magnetic fields across the aperture gives 

H w|f = H j|f + 2H I*U (Over aperture) (1) 

where H- is the incident magnetic field and the subscript t indicates tangential component to the aper- 
ture. Equation (1) is solved for the equivalent magnetic surface currents m, from which the scattering 
properties of the cavity are determined. The electromagnetic fields required to solve equation (1) can be 
obtained by using the procedure described in section 2.1. 

2.1. Electromagnetic Field in Interior Part 

The transverse electric and magnetic fields in the cylindrical cavity may be obtained by the superpo- 
sition of incident and reflected fields as 

P 

E h, = X [ a p exp (-ftp*) + b P exp Uy p z)]e p (x,y) (2) 

p=o 

p 

H W = X [ a p exp - b P ex P UY p z)]Y p h p (x,y) (3) 

P=0 


where a p and b p are amplitudes of the forward and backward traveling waves, y p is the propagation 
constant of pth mode, Y p is the waveguide modal admittance, and e p (x,y) and h^(jt,y) are the 
waveguide vector modal functions as defined in references 10 and 11. The analysis in this paper is lim- 
ited to a rectangular cavity cross section in which fields can be expressed in suitable modal expansion 
functions. Subjecting the transverse electric field at z = -L to zero allows equations (2) and (3) to be 
written as 


P 

e w = x ~ 2 ^ a p exp (ttp L ^ sin e p( x >^ w 

P = 0 
p 

H h; = X la p exp COS Yp{z + L ) Y p h p (x,y) (5) 
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Because at z = 0, 


E,| z=0 x (-2) = -m(*,y) 
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equations (4) and (5) may 


be written as 
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Y p h p (x,y) (7) 


The aperture is now divided into triangular subdomains as shown in figure 2. The triangular discretiza- 
tion shown in figure 2 is arbitrary. The final results will depend upon the average size of the triangular 
elements. If the average size of subdomain triangles is small, we expect the results to be more accurate. 
The triangular discretization shown in figure 2 is obtained by using COSMOS/M software. The mag- 
netic surface current density over the aperture may be expressed in terms of triangular basis functions 
(ref. 13) as 


N 

m(*',y,0) = X T n »„(*•') 
n = 1 


( 8 ) 


where T n is the amplitude of electric current normal to the nth edge, B n (r') is the vector basis func- 
tion associated with the nth edge, and N is the number of nonboundary edges on the aperture surface. 
The expressions for the basis function are given by (ref. 13) 
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where all parameters in equations (9) are defined and explained in detail in reference 13. 

By using equation (8), the transverse electric and magnetic fields from equations (6) and (7) inside 
the cavity may be written as 
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2.2. Electromagnetic Field in Exterior Part 

In the exterior part, the electromagnetic field consists of (1) the field due to the incident wave and 
(2) the field scattered by m(x,y) backed by an infinite ground plane at z = 0. 
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2.2.1. Incident field. With time variation exp (j(Ot) assumed, the incident electric field may be 
written as 




exp (-;'k. • r) 


( 12 ) 


where 


k/ = -*o(* sin 0, cos + y sin 9, sin + z cos 0/) 
r = xx + yy + zz 


E e i = | E ;| cos “o 
E = | E ,-| sin “o 

and &o is the free-space wave number. The angle a Q = 0 corresponds to H-polarization and a Q = 90° 
corresponds to E-polarization. (See fig. 1.) From equation (12), the Cartesian components of the inci- 
dent electric field may be written, respectively, as 


E x = Eq cos 0- cos sin §i 

(13) 

E y = Eq cos 0^ sin + cos 

(14) 

E „ = -Eq sin 0. 

Z : 0; l 

(15) 


The corresponding magnetic field components are obtained through 

H - = i k ‘- x E < (16) 


where ri 0 is the free-space impedance. The incident field with Eq yt 0 and = 0 is called an 
H-polarized wave; with Eq = 0 and E^ ■£ 0, an E-polarized wave. 

2.2.2. Scattered field due to m(jc,y,0). For determining the scattered fields in the exterior part, the 
equivalent magnetic current m(x,y,0) backed by an infinite conducting plane is considered. Consider- 
ing the image of m( ) in the infinite ground plane, the electromagnetic fields for the exterior part may 
be written in terms of the magnetic vector potentials F(m) as 

E„(2m) = -?V x F(m) (17) 

s £ 


H/2m) = -2yco F(m) - 


-^-VV • F(m) 

cojLte 


(18) 


e rf exp [-y^ 0 (r - r')] 

F(m) = s; JJ mt *’- y) — — dS 


( 19 ) 
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where £ and p, are the permittivity and permeability, respectively, of free-space, and r and r' are the 
coordinates of the field and source points, respectively. With equation (8), the scattered electric and 
magnetic fields may be written as 


N 

E s (2m) = X r »7 Vx F ( B „) 

71 = 1 


( 20 ) 
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( 21 ) 


where F(B n ) is computed from equation (18) with B^ replacing m. The total EM fields in the exterior 
region are now obtained by superposition of the scattered and incident fields as 


E t = E,(2m) + £, + £, 


( 22 ) 


H, = H s (2m) + H | + H /? (23) 

where and are the reflected electric and magnetic field parts of the incident fields due to the 
infinite ground plane. In the plane of the infinite ground plane, note that E . + = 0 and 

H ; + H r = 2H,.. 

Substituting equations (11) and (21) into equation (1), the magnetic field integral equation with m 
or T n as an unknown is obtained as 


N P 


I F nX 


~j cos y p L 
sin y p L 


JJ »„(*'•>') 

Aperture 


h p (x\y')dx' d? 


Y p h p (x,y) 


N 


n = 1 


= I r „ - 2 7» F(B„) 


ALvv 

00 (1£ 


F (B„)] + 2H ; | ( 


(24) 


By selecting B /2 , with ri = 1, 2, 3,...,N as a testing function, the application of the method of moments 
to equation (24) yields 
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Equation (24) may be written in matrix form as 

[Y] [r] = [Y] 


(25) 


(26) 
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where the elements of matrix [Y] are given by 


_ cos jL . r 

Y n’n = I sin v j. Y „ J J ' V*>’> * 


sin yL 

P = 0 *P 
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and the elements of [V] are given by 


(27) 


V n , = J J B n ,(x,y) • H ( .| f dx dy (28) 

Equation (26) is used to determine the magnetic current on the aperture. The far field due to m is 
obtained with the procedure in section 2.3. 


2.3. Far Field 

By using the far-field approximation, the scattered electric far fields in the (0,<{)) direction may be 
obtained from equation (20) as (fig. 1) 
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E = -2 V T jkr.F 
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(30) 


where F n and F n are given by 
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47 ir 
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£> 


where r is distance to far field. The copolarized radar cross-section pattern of an aperture backed by a 
cavity in a three-dimensional conducting surface can be obtained from 


a ee 


lim 

f > CO 


47ir 



for H-polarized incidence and 
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(32) 


(33) 
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for E-polarized incidence. The cross polarized radar cross-section pattern is obtained by using 


a <t>0“ 



lim 

47 ir 
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3. Numerical Results 

The RCS of an arbitrarily shaped aperture backed by a cylindrical cavity recessed in an infinite 
ground plane is determined by solving the matrix equation (26) for [Fj. The cylindrical cavity may be 
of rectangular, circular, or coaxial cross section. In fact, the method is applicable to any regularly 
shaped cavity in which EM fields can be expressed in the appropriate modal functions. The matrix ele- 
ments required for the solution of equation (26) are determined from equations (27) and (28) and 
numerical integration over the triangles. For numerical integration over the triangles, a 13-point Gauss- 
ian quadrature formula (ref. 14) is used. Numerical convergence of matrix equation (26) depends upon 
the choices of values of N and P. For numerical convergence, sufficiently large values of N and P were 
selected in all the examples considered. For evaluation of the self terms, that is when n = ri, 
closed-form expressions given in reference 15 are used. The magnetic surface current density obtained 
from the solution of equation (26) is then used in equations (29) and (30). 

A rectangular aperture with x and y equal to 0.7 X Q and 0. 1 Xq , respectively, and backed by a rect- 
angular cavity of depth 1 .73X Q is selected as a first example to validate the present method of moment 
and modal expansion method. This example is selected because its RCS patterns calculated with other 
numerical techniques are available in reference 9. The rectangular aperture for the present method was 
divided into triangular subdomains. The backscatter RCS pattern of the rectangular aperture is com- 
puted as a function of the number of triangular subdomains (i.e., N) and number of cavity modes (i.e., 
P). The results of these computations are presented in figures 3 and 4. The conclusion may be that for 
convergence of numerical results, the average area of the triangular subdomain must be less than 
0.002 X 2 and P > 20 . The cavity modes considered were selected in the ascending order of their cutoff 
frequencies. For the numerical examples considered in this paper, the authors did not observe any spe- 
cific pattern of cavity modes that may enhance the convergence. Figure 5 shows backscatter RCS pat- 
terns of the rectangular aperture computed with the present method and earlier published results of 
reference 9. There is good agreement between the two results. Additional backscatter RCS patterns for a 
long and narrow cavity (a = 2.5 Xq, b = 0.25 X Q ) with depth equal to c = 0.25 Xq computed with the 
present method are shown in figure 6 along with earlier published results of reference 9. 

Figure 7 shows a circular aperture backed by a rectangular cavity recessed in an infinite ground 
plane and excited by an EM plane wave. The RCS pattern of the circular aperture is calculated by the 
present method and presented in figure 8. 

4. Concluding Remarks 

The electromagnetic scattering from a rectangular aperture backed by a rectangular cavity recessed 
in an infinite ground plane has been analyzed by the integral equation approach. The electromagnetic 
fields outside the cavity are obtained in terms of the half-space Green’s function and the magnetic sur- 
face current density flowing over the aperture. The electromagnetic fields inside the cavity are obtained 
in terms of cavity modal functions. The integral equation with the magnetic surface current density as an 
unknown variable is derived by enforcing the continuity of the electric and magnetic fields over the 
aperture. The magnetic current density is determined by solving the integral equation in conjunction 
with the method of moments. The backscatter patterns of rectangular and circular apertures backed by 
rectangular cavities are determined from knowledge of the magnetic current. 
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The numerical results for rectangular apertures backed by rectangular cavities are compared with 
earlier published results. There is good agreement between the results obtained by the present method 
and earlier published data. The radar cross section of a circular aperture backed by a rectangular cavity 
is also presented. Even though numerical results on rectangular and circular apertures backed by rectan- 
gular cavities are presented in this report, the present method can be extended to analyze an arbitrarily 
shaped aperture backed by circular or coaxial cavities. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
December 11, 1996 
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(a) Geometry. 



(b) Equivalent exterior problem. (c) Equivalent interior problem. 

Figure 1. Rectangular aperture backed by a rectangular cavity recessed in an infinite ground plane. 
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Figure 2. Geometry of arbitrary aperture discretized into triangular subdomains. 



= 4>|» de g 

(a) E-polarization. 

Figure 3. Backscatter RCS patterns for rectangular aperture backed by rectangular cavity 
(a = OJX 0 , b = 0.1 X Q , c = 1.73A, 0 ) recessed in an infinite ground plane (0 = 0. = 40.0°) with 
N = 58 and P = 10, 20, 40, and 60. 
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(j) = <j>., deg 

(b) Cross polarization. 
Figure 3. Continued. 
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IX 2 , dB 





Figure 4. Backscatter RCS patterns for rectangular aperture backed by rectangular cavity 
( a = 0.7 Xq, b = 0.1 Xq, c = 1.73A,q) recessed in an infinite ground plane (0 = 0 ; . = 40.0°) with 
P = 40 and N = 22, 48, 58, and 108. 




10 ,- 



<j) = <j)., deg 

Figure 5. Backscatter RCS patterns for rectangular aperture backed by rectangular cavity 
( a = 0.7A-Q, b - 0.1 A,q, c - 1.73Xq) recessed in an infinite ground plane (0 = 0. = 40.0°) with 
P = 40 and N = 58 and results from reference 9. 
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(a) <j> = <j>. = 0.0°. 

Figure 6. Backscatter RCS patterns for rectangular aperture backed by rectangular cavity 
(i a = 2.5 A,q, b = 0.25 A,q, c = 0.25 A, 0 ) recessed in an infinite ground plane with P = 40 and N = 108 
and results from reference 9. 
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(b) <J) = <j). = 90.0°. 
Figure 6. Concluded. 
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